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ABSTRACT 

We have studied the stability of low degree g-modes in uniformly rotating B-type stars, 
taking into account the effects of the Coriolis force and the rotational deformation. 
From an analysis treating rotation frequency as a small parameter it is found that slow 
rotation tends to destabilize high radial-order retrograde g-modes, although the effect 
is very small or absent for relatively low order modes. Calculating eigenfrequencies 
at selected rotation rates, we find, on the other hand, that rapid rotation tends to 
stabilize retrograde (/-modes. The stabilizing effect appears stronger for less massive B- 
type stars having low effective temperatures. If we change rotation rate continuously, 
the frequency of a g-mode belonging to (£, m) crosses frequencies of other g-modes 
belonging to (£',m). If the parity of the two encountering modes are the same, they 
interact each other and the stability (i.e., imaginary part of eigenfrequency) of each 
mode is modified. Using an asymptotic method we discuss the property of such mode 
crossings and couplings. For rapidly rotating stars mode couplings are important for 
the stability of low degree g-modes. In particular, we find that the stabilization of 
retrograde g-modes in rapidly rotating stars is due to many strong mode couplings, 
while prograde sectoral modes are exceptionally immune to the damping effects from 
the mode couplings. 
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1 INTRODUCTION 

The kappa-mechanism associated 
ity bump at T ~ 2 x 10 5 K 
modes in intermediate-mass (7 
main-sequence stars 



with the iron opac- 
excites high-order 
> M/Mp > 3 
dGautschv fc Said [l99 



19931 ). 



These B- 



iDziembowski. Moskalik. Pamvatnvkh 
type g-mode pulsators with typical periods of the order 
of days a r e call ed slowly pulsating B (SPB) stars, which 
IWaelkensI i|l99ll ) first recognize d as an independent class 
of variables (see |Pe Cat! 120071 . for a review). It is gen- 
erally considered that majority of SPB stars are slow 
rotators with projected equatorial rotation velocities of 
Usini < 100 km s _1 , although several SPB stars are known 
to have Usini exceeding 200 km s" 1 . However, the criterion 
for slow rotation may not be appropriate when we look into 
the effect of rotation on g-modes of SPB stars. For example, 
an equatorial rotation velocity of 50 km s _1 corresponds 
to an angular velocity f2 — 0.1a/ GM/R 3 if we adopt mass 
and radius typical for a SPB star; M = 4M Q and radius, 
R « 2ARq, where G is the gravitational constant. This 
indicates that the centrifugal acceleration, RQ. 2 is much 



* E-mail: aprilia@astr.tohoku.ac.jp 
f E-mail: lee@astr.tohoku.ac.jp 
| E-mail: saio@astr.tohoku.ac.jp 



smaller than the surface gravity, GM/R 2 , causing negligible 
deformation of the structure of the star. However, the 
50km s _1 velocity corresponds to a rotation period of about 
2.4days, comparable to the typical pulsation periods of 
SPB stars; in other words, the ratio of rotation to pulsation 
frequency is comparable to unity. Since the Coriolis force 
affects stellar oscillations significantly when 2Q/u > 1 with 
uj bei ng the oscillation frequency in the co-rotating frame 
(e.g., lUnno et all 1 19891 ). g-modes in most SPB stars might 
be affected significantly by the Coriolis force, even though 
rotational deformations are negligibly small. 

In addition, it is known that g-modes are 
present in many rapidly rotating Be stars (e.g., 
iRivinius. Baade. fe Stefll 120031 ). This is understandable 
because many Be stars are located in the SPB instability 
strip on the HR diagram. There is some evidence that 
the excitation of g-modes in Be stars should be affected 
signifi cantly by rotation (e.g.. IWalker et alT 2005: Sai o et all 
120071 ). Therefore, it is important to examine the effect of 
rotation on g-modes in B-type stars. 

Stability analyses of g-modes in rotating B s tars 
have been carried out by several authors (e.g. , iLee 



200ll; iTownsendl 120051; ISavoniid [20051; IWalker et alj 12005. 
Saio et al.ll2007l ; ICameron et aU2008r i.lLeej i|200ll ) calculated 



non-adiabatic g-modes of rotating SPB stars by taking into 
account the effect of the Coriolis force and found that the 
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low degree unstable retrograde g-modes in the absence of ro- 
tation are stabilized by ra pid rotation via mode couplings 
with high i stable o-modes. lSavoniiel (|2005l ) obtained similar 
results on the stabilization of retrograde modes in rapidly ro- 
tating models, by solving a set of two dimensional partial dif- 
ferential equations describing small amplitude non-adiabatic 
oscillations of uniformly rotating stars. The tendency of ret- 
rograde modes being stabilized is more pron ounced for g- 
modes in nearly critically rotating Be star s (| Walker et al.l 
120051 ; ISaio et al.ll2007l ; ICameron et al1l2008h . 

On the other hand, such stabilization of retrograde 
modes does not occur if the traditional approximation is em- 
ployed (|Townsendll2005l ). The horizontal component of the 
angular velocity of rotation (Q. sin 9 with 6 being co-latitude) 
is neglected in the traditional approximation, which makes 
the set of the governing equations similar to the one with- 
out rotation. Only difference is that 1(1 + 1) in the non- 
rotating case is replaced with \k m (a function of 20. /lj) in 
the traditional approximation, where k = 0, 2, 4, . . . for even 
modes, while k = 1, 3, 5, . . . for odd modes, in parallel to 
I = |m| + k in the non-rotating case. The traditional ap- 
proximation works well for g-modes except when two modes 
associated with the same m and parity (even or odd) but 
different Afe m s have similar frequencies; in the traditional ap- 
proximation the two mo des are independen t, while the two 
modes couple in reality ijLee fc Said 1 19891 ). Therefore, the 
difference in the stability result based on the traditional ap- 
proximation from the result without it can be explained by 
the absence of mode couplings under the traditional approx- 
imation. This indicates that mode couplings are important 
in the stability of g-modes in rotating stars. In this paper, we 
discuss the effect of the couplings on the stability of g-modes 
in detail. 



2 METHOD OF CALCULATION 



a 



We u se the method of calculation given in iLee fc Baraffd 
i|l995l ) to study the stability of g-modes in uniformly rotat- 
ing stars by taking into account the effects of the Coriolis 
force and the centrifugal force. We employ a coordinate sys- 
tem (a, 9, (j)) for a rotationally deformed star, where the co- 
ordinate a is regarded as the mean distance of equi-potential 
surface from the stellar centre. It is related to spherical polar 
coordinates (r, 8, (j>) as 



r = a[l + e(a,0)], 

where e, a rotational deformation, can be written as 
e = a(a) + /3(a)P 2 {cos 9) 



(1) 



(2) 

with the second Legendre polynomial P2(cos 9) — (3 cos 2 9 — 
1) /2. The term a represents the spherical expansion and 
P the deformation of the star due to rotation. Assuming 
that e is proportional to Q. 2 , we calculate the functions a(a) 
and P(a) by applying the Chandrasekhar-Milne expansion 
to the hydrostatic and Poisson equations for the star (see 
ILee fc Baraffelll995l . for details). 

We express the angular dependence of small amplitude 
oscillations of the star using finite series expansion in terms 
of spherical harmonic functions Y ; m (#,0). For a given az- 
imuthal wavenumber m, the displacement vector £(a, 9, <j>, t) 
is given by 



3=1 

jmax 

E 



3=1 



3=1 



H h {a)\ 



SrYiT - Tu (o) ^F;' 



(3) 



(4) 



(5) 



and the Eulerian perturbation of the pressure, p' , for exam- 
ple, is given by 



(a,W) = ^P;>)^7(MK 

3=1 



(6) 



where lj = a + mO is the oscillation frequency in the co- 
rotating frame of the star with a being the oscillation fre- 
quency in an inertial frame, and lj = jmj + 2(j — 1) and 
l'j = + 1 for even modes, and lj = \m\+2j — 1 and lj — — 1 
for odd modes, where j — 1, 2, ■ • ■ , j max , and j max is the 
length of the expansions. Note that the angular pattern of p' 
for even (odd) modes is symmetric (antisymmetric) with re- 
spect to the equator. Substituting these expansions into the 
linearized basic equations, we obtain a finite set of coupled 
first order linear ordinary differential equations for the ex- 
pansion coefficients (a) and p[. (a), and for given values of 
the parameters m and fl, we solve the set of the differential 
equations as an eigenvalue problem of uj by imposing appro- 
priate boundary conditions at the centre and surface of the 
star. The set of the linear differe ntial equations as wel l as the 
boundary conditions is given in ILee fc Baraffel(|l995l ). Since 
we assume that the perturbed quantities are proportional to 
e i(m<t>+uit) ^ moc [ es nav i n g negative Wi are pulsationally unsta- 
ble (or excited), where ui\ = lm(uj) and ojr = Re(o;) denote 
respectively the imaginary part and the real part of lj. Note 
that when cJn > 0, modes associated with negative (pos- 
itive) azimuthal wavenumber m are prograde (retrograde) 
modes seen in the co-rotating frame of the star. 

To save computing time, we employ the Cowling ap- 
proximation, neglecting the Eulerian perturbation of the 
gravitational potential. We ignore the effect of rotational 
spherical expansion of the equilibrium structure on the os- 
cillation, that is, we set a = in equation because it is 
just a change in mean radius not important for our qualita- 
tive study. 

To determine an adequate length of the series expan- 
sions, j m ax, we have compared the results of sample calcula- 
tions obtained with j ma x = 6, 8, 10, and 12. We have found 
that jmax = 10 is large enough for un stable ff- modes even 
for the case of rapid rotation Q, = 0.4 GM/E?. Therefore, 
we have adopted jmax = 10 in calculating eigenfrequencies 
and eigenfunctions of g-modes presented in this paper. 



2.1 Mode identification and notation 

As equations <[3j — <[0J indicate, the angular variation of a non- 
radial pulsation mode in a rotating star cannot be repre- 
sented by a single spherical harmonic but need to be ex- 
panded into many harmonics for a given azimuthal order 
m. The angular variation of a mode depends on the ratio 
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Figure 1. Eigenvalue of Laplace's tidal equation A^ m versus 
— sgn(m)v for \m\ = 1, where sgn(ra) = m/\m\ and v = 2Q/u). 
Solid lines are for even modes and dashed lines for odd modes. 



Table 1 . Physical parameters of the main sequence models with 
X = 0.7 and Z = 0.02 



Model 


l0g 4 


M = 4M 
logT cff 


X c 


M 


= 5M 
logT cff 


x c 


A 


2.371 


4.165 


0.700 


2.727 


4.226 


0.700 


B 


2.464 


4.144 


0.476 


2.825 


4.208 


0.481 


C 


2.519 


4.120 


0.308 


2.886 


4.185 


0.311 


D 


2.565 


4.083 


0.105 


2.939 


4.149 


0.105 



and lj = \m\ + 2(j — 1) and l'j = l 3 ■ + 1 for even modes, and 
I j r = jmj + 2j — 1 and = lj ■ — 1 for odd modes. 

If we find fj ^0.5 with a particular value of j for a cal- 
culated even mode with a given m, we identify it as a mode 
associated with A^ m with k = 2(j — 1) (if the mode is an 
odd mode k = 2j — 1). For the case of rapid rotation, how- 
ever, we sometimes need to use a somewhat relaxed criterion 
given by fj a with a < 0.5, particularly for modes with 
j > 1. Also, we denote the mode as (£,m) with I = \m\ + k 
extending the familiar notation for nonradial pulsations of 
a non-rotating star. Note that prograde sectoral modes cor- 
respond to £ = — m. 



of rotation frequency to the oscillation frequency. To dis- 
tinguish the angular dependence property of each mode, 
it is useful to refer to the mode properties under the tra- 
ditional approximation, in which the term — fi sin Oeg in 
ft = £1 cos 6e r — Q sin 6eg is neglected. Many properties of 
low frequency oscillations of uniformly rotating stars are 
well explained by using the tr aditional approximation ex- 
cept for mode couplings (e.g., iLee fc Saiol Il987l . 119971 ). In 
this approximation, the angular depend ence of the oscilla- 
tions are given by Hough function (e.g., iLindzen fc Holtonl 
1968) which is the eigenfunction, associated with the eigen- 
value Xkm, of Laplace's tidal equation. For a given azimuthal 
wavenumber m, Afc m depends on the ratio v = 2Q,/uj and 
tends to 1(1 + 1) with I = \m\ + k as v -> for k ^ 0. Fig. [1] 
shows Xkm as a function of v for k — — 3 with \m\ = 1. The 
quantity \/Xk m represents a kind of surface wave number. 
Except for the prograde sectoral modes (A: = 0), Xkm in- 
creases as the parameter v increases. The prograde sectoral 
modes (associated with Xo m (m < 0)) are special modes in 
rapidly rotating stars whose surface wavenumber is lower 
than the value at = 0, and hardly changes with £1. We 
emphasize here that although we use the modal properties 
under the traditional approximation as guidelines to under- 
stand the behavior of g-modes as function of f2, we do not 
use the traditional approximation in our numerical analyses. 

To identify calculated modes, we make use of the frac- 
tional kinetic energy of oscillation mode defined as 



(7) 



where 



+l'j (fj + 1) 





2' 






3 





in pa 4 da, 



(8) 



2.2 Unperturbed models 

For unperturbed models of the stability analyses, we have 
calculated main-sequence evolution models without rota- 
tion, usi ng a standard stellar evolut i on code using the OPAL 
opaci ty (|lglesias. Rogers fc Wilson! 1 19921 ; Egj esias fc Rogers! 
1996). We have adopted masses of 4Mq and 5Mq with ini- 
tial chemical compositions of (X, Z) = (0.70, 0.02) as typical 
models for SPB stars. Table [T] shows parameters of selected 
models. 



3 G-MODES EXCITED AT SELECTED 
ROTATION FREQUENCIES 

Figs. EHS] show growth rates, r\ = — wi/wr, of unstable g- 
modes in the 4M and 5M ZAMS and D models (Tab.[TJ 
at three selected rotation frequencies; Q. — 0.01, 0.2, and 0.4, 
where Q = Q/y/GM/R 3 . The horizontal axes — sgn(m)wR 
measure the oscillation frequency in the co-rotating frame 
of the prograde (retrograde) modes in positive (negative) 
direction, where lj — ujjyl GM/R 3 . Cyan and red dots are 
for (£, |m|) = (1, 1) and (2, 1), respectively, while crosses and 
black dots are for (I, |m|) = (2, 2) and (3, 2), respectively. 

Firstly let us look at the slowest rotation cases (left pan- 
els). When rotation is very slow, crosses and red dots are al- 
most superposed on each other, because eigenfrequencies are 
independent of the azimuthal order m in a non-rotating star, 
and the first order rotational correction to the frequency, 
which is proportional to m, is small for slow rotation. The 
growth rates n as a functions of ljb, are almost symmetric 
with respect to the axis of lDr = 0. Many gr-modes are ex- 
cited in a certain frequency range depending on the effective 
temperature of each model; the unstable frequency range is 
lower for a lower effective temperature, which is obvious, 
for example, when a ZAMS case is compared to the case of 
the evolved D model. (We note that the density of g-modes 
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Figure 2. Growth rate r) = — cjj/cjr versus — sgn(m)iI)R for g-modes excited in the 4Mq ZAMS model A, where the cyan, red, and black 
dots stand for the (/-modes with (£, \m\) = (1, 1), (2,1), and (3,2), respectively, and the crosses stand for the (/-modes with (£, \m\) = (2, 2). 
Left-, middle- and right-panels are for rotation frequencies of Q = 0.01, 0.2 and 0.4, respectively. Prograde modes having m < appear 
on the right half of each panel, while retrograde modes (m > 0) appear on the left side. 
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Figure 3. Same as Fig. [2] but for the 4Mq evolved main sequence model D. 
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Figure 4. Same as Fig.[2]but for the 5M Q ZAMS model A. 
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Figure 5. Same as Fig. [2] but for the evolved 5Mq main sequence model D. 



is higher in the latter because the Brunt- Vaisala frequency 
increases as evolution proceeds.) This tendency of stability 
arises from the optimal condition for the kappa-mechanism 
of excitation; the kappa-mechanism is most effective when 
the thermal-timescale at the Fe opacity bump is comp arable 
to pulsation periods (see e.g.. iGautschv fe Saiolll993l ). That 
is, longer-period modes tend to be excited in cooler models 
in which the opacity bump is located in deeper layers having 
longer thermal-times. 

In the 5M Q ZAMS model (Fig. |4j , £ = 1 g-modes are 
not excited at Q — 0.01, because their frequencies are too 
low to satisfy the optimal frequency condition for the kappa- 
mechanism. Retrograde I = 1 g-modes become unstable as 
the rotation frequency increases slightly as seen in the mid- 
dle panel of Fig. 0] This is because the frequencies of ret- 
rograde modes have increased to enter into the frequency 
range where the kappa-mechanism works strongly enough to 
excite modes. Since the rotational effect on the frequency in 
a very slowly rotating star is approximately given as mCiQ 
in the co-rotating frame with a positive constant C\ (see 
the next section), frequencies of retrograde modes increase 
as the rotation rate increases. 

It is seen in the panels for f2 = 0.2 and 0.4 in Figs. [2]- 
[S] that the prograde-retrograde symmetry is broken as the 
rotation frequency increases; the maximum growth rates for 
the retrograde modes become lower than those of the pro- 
grade modes, and the oscillation frequencies of the unsta- 
ble retrograde modes increase as the rotation frequency in- 
creases, while the effect is not appreciable for the prograde 
modes. We also note that the number of unstable g-modes 
is largely reduced, particularly for the retrograde g-modes. 
This is due to selective damping caused by mode couplings, 
which will be discussed in iJS] 

Fig. [5] shows unstable g-modes of a 5Mq evolved main- 
sequence model, which has an effective temperature similar 
to that of the 4M Q ZAMS model (see Table 1). Comparing 
Fig. [5] with Fig. fj for the the 4M ZAMS model, we find 
that although the frequency spectra for the former model 
are much denser than those of the latter, rapid rotations 
stabilize retrograde modes similarly for both models. This 
may suggest that the stabilizing effect on retrograde g-modes 
mainly depends on the effective temperature, while the dif- 



ference in the internal structure due to mass difference has 
only a minor effect. 

Figs. [2H5] clearly show that in rapidly rotating stars 
retrograde modes tend to be stabilized more strongly than 
prograde modes. In particular, for the 4M evolved model 
shown in Fig. only prograde (mostly sectoral I = 2) modes 
are excited. The retrograde-prograde mode asymmetry is 
less pronounced in relatively hotter models. This is (at least) 
partly because unstable frequency ranges in hotter mod- 
els are higher (due to the optimal condition for the kappa- 
mechanism at the Fe opacity bump) so that the ratios 20. /ui 
(which govern the strength of the Coriolis force effect) for 
the excited modes are generally smaller, and hence the ro- 
tation effects are weaker in hotter models. 



4 SLOW ROTATION 

Let us next discuss the effect of slow rotat ion on th e sta- 
bility of g-modes in rotating B stars (e.g., iLeej [200l|). For 
oscillations in a slowly rotating star, we can treat the ro- 
tation frequency SI as a small parameter, and express the 
oscillation frequency uj of a mode as 

Q (fi) = Cjq + mCiQ + C 2 S1 2 , (9) 

where cjo is the oscillation frequency of the mode in the ab- 
sence of rotation. The coefficient C\ represents the first order 
rotational effect due to the Coriolis force and the coefficient 
C2 the second order effects that come from the centrifugal 
force and the Coriolis force. Since eigenfrequencies tj are 
complex numbers in nonadiabatic analyses, C\ and Ci are 
also complex numbers. The real part of C\, Cm, is approx- 
imately equal to the adiabatic expression 

ClR ~ / *[2WU + l]&pr 2 rfr 
1R ~ J R [(WU) 2 + e(£ + mLpr 2 dr 

(e.g.. lUnno et al1ll989l ). where subscript 'a' indicates adia- 
batic eigenfunctions. 

We note that negative (positive) Cn = Im(Ci) means 
that the Coriolis force due to a slow rotation has destabi- 
lizing (stabilizing) effect on retrograde modes with m > 
and stabilizing (destabilizing) effect on prograde modes with 
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m < 0, and that negative (positive) C21 means the desta- 
bilizing (stabilizing) effect on both prograde and retrograde 
modes. In this paper, we have obtained the complex coeffi- 
cients G\ and C2 for a given m, by computing the eigenfre- 
quencies to of a mode at three different rotation frequencies, 
e.g., H = 0, 10 -3 and — 10~ 3 , and substituting these val- 
ues in equation ([9]). We confirmed that Ci thus computed is 
in good agreement with that obtained with the method by 
ICarroll fc Hansen! (|l982T ). 

In Fig. [51 we plot Cl and C2 versus ijr for £ = \m\ — 1 g- 
modes of 4M© main sequence models, of which physical pa- 
rameters are tabulated in Table [JJ For the ZA MS model C l 
and C2 are essentially the same as the result of I Led (|200lT ). 
The coefficients Cm. for the evolved models are systemat- 
ically larger than that for the ZAMS model. This can be 
understood as follows. In the interior of an evolved model 
the Brunt- Vaisala frequency N is higher than in the ZAMS 
model. An asym ptotic theory of nonradial pulsations (e.g., 
IShibahashi|[l979l ) gives the ratio of radial to horizontal dis- 
placement as ~ yJl(J+V)\uj /N\, which indicates the 
ratio |£r/£fc| is smaller in the evolved models. The smaller 
ratio, in turn, yields larger Cm as seen from equation (110[) . 

In addition, Cm for an evolved model has many reg- 
ularly spaced peaks. The peaky features are also found in 
the imaginary parts of C\ and C2 for the evolved models. 
We have also calculated C\ and C2 for WMq main sequence 
models (not shown) to find quite similar behavior of the coef- 
ficients as functions of uir. The peaks in Ci are caused by g- 
mode eigenfunctions being trapped in the /i-gradient zone (/x 
stands for mean molecular weight) above the convective-core 
boundary, in which the Brunt- Vaisala frequency is much 
higher than the adjacent layers. As discussed above, this 
yields larger Cm for the trapped modes. 

The regularity of the peaks is also understood from the 
asymptotic theory, which gives the frequency of a (/-mode 
trapped in the /x-gradient zone, as 



N 



dr. 



(11) 



where Ar and n M denote respectively the width of the \i- 
gradient zone and the number of radial nodes within the 
zone. Equation indicates that the frequency separation 
between two consecutive trapped (/-modes is proportional 
to [n M (n M + 1)] _1 , which explains why peaks are reqularly 
spaced and the separations of peaks decrease with increasing 
radial orders (or decreasing frequencies). 

The imaginary parts of Ci and C2 decrease rapidly as 
ujb, decreases in the range uj-r < 0.3. Since the imaginary part 
of Ci is negative for high radial-order (/-modes, slow rotation 
has destabilizing (stabilizing) effect on high radial-order ret- 
rograde (prograde) (/-modes. Considering that the imaginary 
part of C2 is also negative for the high radial-order (/-modes, 
the destabilizing effect of slow rotation on the (/-modes works 
more strongly for the high radial-order retrograde modes. 
For relatively low radial-order modes both Cn and C21 are 
very small, indicating a slow rotation affects little the sta- 
bility of these modes. The steep decreases of the imaginary 
parts of Ci and C2 with decreasing g-mode frequencies can 
be understood from the relation between the period of a 
high order g-mode and the optimal period for the kappa- 
mechanism excitation that is roughly the thermal timescale 



at the Fe opacity bump. The period of a very high-order 
g-mode that has negative Cn and C21 is much longer than 
the optimal period for the kappa-mechanism. If the mode 
is a retrograde mode, both the first order and the second 
order effects of rotation increase the frequency (or decrease 
the period) as eq.([5| indicates. This effect makes the kappa- 
mechanism to work stronger to the mode, and appears as 
negative Cn and C21; in other words, the first and the sec- 
ond order rotation effects tend to destabilize very high-order 
retrograde g-modes. On the other hand, for a prograde mode 
the first order effect increases the period, while the second 
order effect decreases the period. Therefore, the first order 
rotation effect tends to stabilize a prograde mode, while the 
second order effect tends to destabilize it. 



5 MODE COUPLINGS 
5.1 Numerical Results 

Fig. [7] shows the complex frequencies uj of selected even g- 
modes of (£, \m\) — (1, 1) (associated with Ao±i) (red lines) 
and of (£, \m\) — (3, 1) (associated with A2±i) (black lines) as 
functions of O for the 4Mq ZAMS model, where the radial 
orders of the (/-modes range from 21 to 26 for the £ = 3 
modes and from 9 to 11 for the £ — 1 modes. The left panels, 
(a) and (b), are for prograde (m = —1) modes, while the 
right panels, (c) and (d), are for retrograde (m = +1) modes. 

The variation of pulsation frequency ljr of a g-mode 
as a function of reflects the variation of Xkm as a func- 
tion of v. For a high order g-mode, pulsation frequency is 
approximately written as 




(12) 



where n is the radial order, N = N/^fGM/R 3 normal- 
ized Brunt-Vaisal a frequ ency, and x fractional radius r/R 
l|Lee fc Said 1 19871. Il989h . The equation indicates that the 
frequency of a prograde sectoral mode (£ = — m = 1) 
slightly decreases but stays nearly constant as the rotation 
frequency increases just as A0-1 behaves as a function of 
v. The frequencies of the other modes in Fig. [7] increase at 
different rates depending on the corresponding Xk m s. There- 
fore, the oscillation frequency of a (/-mode associated with 
a given Afc m crosses with another g-mode associated with 
Afe' m (k' 7^ fc) as seen in Fig. [JJ 

For adiabatic modes having real frequencies, a mode 
crossing between (/-modes leads to an avoided crossing, but 
for non-adiabatic modes having complex frequencies, mode 
crossing is not always an avoided one. More importantly, a 
mode crossing between non-adiabatic (/-modes can change 
the stability of the modes. The imaginary parts of eigenfre- 
quencies as functions of Q. in Fig. [JJ show numerous narrow 
peaks large and small. We consider that these peaks are 
caused by mode couplings. In panel (d) of Fig. [7J we recog- 
nize two pairs of a small peak and dip around f2 ~ 0.06 and 
~ 0.07, while panel (c) shows that a black line crosses a red 
line at each of the same H, indicating these peaks and dips 
in ui\ are caused by mode crossings. Large peaks in u\ can 
be interpreted as couplings with strongly damped modes as- 
sociated with higher radial orders and higher A fc / m having 
fc' = 2, 4,---. 
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Figure 6. Complex coefficients C\ and C2 versus cjfi for I = \m\ = 1 g-modes of the 4Mq main sequence models, where the black, red, 
and blue lines indicate, respectively, the models A (ZAMS), B, and C in Table ITl 



Since |wi| of the £= 1 (Ao±i) g-modes is much smaller 
than that of the I = 3 (A2±i) (/-modes, zoomed diagrams of 
the former case are given in Fig. [5] For the I = 1 retrograde 
modes (panel (b)) Qi behaves similarly to I ^ 3 cases with 
numerous high peaks. This is reasonable because the fre- 
quency ub. of a retrograde 1=1 mode increase with Q. sim- 
ilarly to the modes of £ ^ 3 so that similar mode couplings 
would occur. On the other hand, the stability of prograde 
£ = 1 (A0-1) modes behaves differently; for f2 > 0.2, sharp 
peaks are replaced with broad bumps which become weak 
as f2 increases, and those unstable modes stay unstable even 
at large Q. This property is common for prograde sectoral 
modes and consistent with the results in Figs. [5] and [3] which 
show prograde sectoral modes being excited even at a rapid 
rotation. 

Let us now look into two crossings which occur between 
modes of (£,m) = (1, 1) and (3, 1) (or between A01 and A21) 



in detail. Fig. [5] shows a zoomed diagram around a cross- 
ing at Cl ~ 0.06 which is recognized as connected small 
peak and dip in Fig. [7] (panel d). The crossing occurs be- 
tween the modes gg of (£,m) = (1,1) (red line) and (722 
of (£, m) = (3, 1) (black line). As this figure shows, the real 
parts of the eigenfrequencies make an avoided crossing, while 
the imaginary parts make a true crossing. At the crossing 
the properties of the modes are interchanged. 

Fig. 1 101 shows another example of mode crossings which 
occurs between t/io-mode of (£,m) = (1,1) (red line) and 
525-mode of (£, m) = (3, 1) (black line) at fi w 0.026. The 
real parts of the frequencies appear to make a true crossing, 
while the imaginary parts show only a very small dip and 
bump. This crossing may be understand as a crossing with 
an interaction much weaker than the case of Fig. [9] The weak 
interaction reduces the distance at the closest encounter of 
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wrs to zero and causes only very small effects on the stability 
of the modes. 



5.2 Asymptotic Analysis 

5.2.1 Coupling coefficient 

It is helpful to discuss the mode crossing phenomena by 
using an asymptotic method based on the traditional ap- 
proximation. Under the traditional approximation, the os- 
cillation modes of a rotating star are separated into indepen- 
dent modes associated with different Afe m s. The traditional 
approximation works well for p-modes except when the fre- 
quencies of two modes with the same m and the same (even 
or odd) parity become very close to each other. Under the 
traditional approximation these two modes are independent, 
but in reality the term — QsinOeg, neglected in the approx- 
imation, brings about mode coupling between them. 

Here, we extend the asymptotic analysis for adiabatic 



oscillations developed by I Lee fc Said (| 19891 ) to the case of 
non-adiabatic oscillations. Note that no effects of the cen- 
trifugal force are included in the following asymptotic treat- 
ment , and only the effe cts of the Coriolis force are consid- 
ered. iLee fc Said |l989j) have derived a dispersion relation 
which should be satisfied by adiabatic g-modes associated 



tan tan ^ = 



where 



*, ; = - / kjdx + -f 7T, 



(13) 



(14) 



x c denotes the outer boundary of the convective core, n e is 
the effective polytropic index at the surface, and 



k 3 = (15) 
x ui 

The right hand side of equation (|13|) . ejk, represents the 
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Figure 8. Zoomed diagrams for u>\ of even I = \m\ = 1 (/-modes of the 4Mq ZAMS model as functions of f2 for prograde modes in panel 
(a) and for retrograde modes in panel (b), where the black, red, and blue curves denote the gg, gig, and cjn-modes, respectively. 



effect of coupling between the two g-modes. We call it "cou- 
pling coefficient". Note that if the coupling were absent (i.e., 
tjk = 0), equation (|13|) would be reduced to tan^j = and 
tan ^ = 0, which give essentially the same expressions for 
Qj and u)k as equation (If 2[) . 

The coupling coefficient tjk arises from a deviation from 
the traditional approximation aroun d a crossing b etwee n g- 
modes associated with A jm and Xkm. [Lee fe Said (|f 989T ) de- 
rived the expression 



where 



[£f( 


\jrn 


1/4 









Xjf, 



Gj k cos x 3 sinx 

1/4 



(16) 



kjdx + —tv 



(17) 



and Gjk is a matrix element denned in lLee fc Said Hi IS!)), 
which consists of terms conflicting with the traditional ap- 
proximation and brings about coupling between the two 
modes. The matrix element Gjk is proportional to v so that 
€jk — > as v — > 0. This corresponds to the fact that modes 
with different Is are independent in a non-rotating star. 

From equation ()16|) . we can calculate tjk as a function 
of v (or Co) for given m, Q, and (j, k). However, only discrete 
values of v satisfy the dispersion relation in equation (|f 3[) . 
We obtained frequencies satisfying the dispersion relation by 
solving equation (|f 3[l numerically with Newton's method in 
the frequency range between £D = 0.1 and Co — 1 for the 4Mq 
ZAMS model. The values of tjk at the discrete frequencies 
thus obtained are shown in Fig. [TT] for (j, k) = (0, 2), and in 
Fig. [l2]for (j,k) = (0,4) and (2,4). Each point corresponds 
to a (/-mode belonging to Xj, m or Xk, m - Since Aj, m (and/or 
Xk,m) increases rapidly with v (except for the case of Aom 
with m < 0; Fig. [!}, a large number of g-mode frequencies 
enter into the range 0.1 ^ ui ^ 1 (cf. eg. I12|) . 



Fig. If f I shows the coupling coefficient tjk between even 
gr-modes associated with Ao m (£ = \m\) and A2 m (£ ~ \m\+2) 
for \m\ = 1 (panel a) and |m| = 2 (panel b) as functions of v. 
The coupling coefficients increase rapidly as the parameter v 
increases, which is reasonable because without rotation two 
modes with different is axe independent from each other. 
For prograde (m < 0) modes, the coefficients level off for 
v > 2 so that the coupling coefficients of prograde modes are 
appreciably smaller than those of retrograde modes when v 
is sufficiently large. 

Figure [T2l shows coupling coefficients between g- modes 
associated with Ao±i (£ = 1) and A4±i (I — 5) in panel (a), 
and between those associated with A2m {I = 3) and A4 m 
(£ — 5) in panel (b). Obviously, the latter is much larger 
than the former at a given v. 

Fig. II lb . shows that coupling coefficients are larger if 
H = 0.4 (red dots) is used rather than 0.2 (black dots), 
indicating coupling is stronger in a rapid rotator for a given 
value of v. 

We also found that the coupling coefficients for the 5Mq 
ZAMS model (not shown) are similar to those of 4M Q ZAMS 
model shown in the figures. 



5.2.2 Quasi- adiabatic extension 

In order to study mode couplings between nonadiabatic g- 
modes with the dispersion relation, it is necessary to ex- 
tend it by including non-adiabatic effects. We modify the 
wavenumber kj in the quasi-adiabatic approximation as 




1 + ift 



St (V ad - V) 



3 S T (V ad - V) + V„ 
l|Leelll985l ), where the oscillation frequency 



(18) 



u is now re- 



garded as a complex number, and 
TV 2 1 A — 3 



A, 



CO 2 UJC2W ' 



01 



\Kr A pTc v j GM 



(19) 
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Figure 9. Mode crossing between I = m = 1 gg mode (red line) 
and I = 3,m = 1 322 mode (black line). The real part makes an 
avoided crossing, where the mode property is exchanged at the 
closest approach. 



S T = -{dhxp/dhxT) p , V ad = (ainT/91np) ad , V = 
dlnT/dlnp, = dlnp/dlnp, V = — dlnp/dln r, and 
L r is the radiative luminosity and c p is the specific heat 
at constant pressure. Here, for simplicity we have consid- 
ered only radiative damping as the non-adiabatic effect. We 
note that if we define the local thermal time scale Tth as 
T~th = 4ivr 2 H p pTc p / L r with H p — — dr/d\np being the pres- 
sure scale-height, we have UJC2/V = 2nrth/P with P — 2tt/uj 
being the oscillation period, and the term containing fjj, 
which represents the radiative damping effect, is propor- 
tional to PI Tth- For the quasi-adiabatic approximation to 
be valid, we take account of the non-adiabatic contribution 
only in the region where P/Tth <0.01. This criterion, how- 
ever, is somewhat ambiguous since we ignore destabilizing 
contributions near the surface where P/rth ~ 1- 

Replacing the radial wavenumber in equation (|15[) by 
the complex wavenumber given in equation fTSJ, we solve 
the dispersion relation in equation (|13[) to obtain complex 
frequencies ui as functions of fi, where we use n e — 3.525. 
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Figure 10. Mode crossing between (£, m) = (1,1) 310 (red line) 
and (£, m) = (3, 1) 925 (black line) modes. 



Note that since we include only radiative damping as the 
non-adiabatic contribution, we obtain only stable modes 
having positive U\. 

Figs. H3H15l show three examples of mode crossings cal- 
culated using the dispersion relation in equation (|13|) . Fig. 1131 
corresponds to the mode crossing in Fig. [9] in which the real 
parts of the eigenfrequencies make an avoided crossing and 
the imaginary parts make a true crossing causing a large 
peak and dip. Fig. [14] corresponds to that in Fig. 1101 in 
which the real parts of the eigenfrequencies appear to make 
a true crossing and the imaginary parts show a small dip and 
bump. The coupling coefficient tjk for the case of Fig.[l3]is of 
the order of 10 -4 , while for the case of Fig. fH]it is of the or- 
der of 10 -6 , supporting the previous conjecture that a weak 
interaction leads to such a crossing as seen in Fig. [10] The 
value of Q at which the mode crossing takes place does not 
necessarily agree well between the asymptotic method and 
the pulsation calculation. A possible reason for the disagree- 
ment may be attributable to neglecting the effect of the cen- 
trifugal force in the asymptotic treatment. The derivatives 
duju/dfl for low frequency g-modes are different between the 
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Figure 11. Coupling coefficient ej^ between even <j-modes asso- 
ciated with A() m and A2m for \m\ = 1 in panel (a) and for \m\ = 2 
in panel (b) as a function of v = 2Q/u. Here, we have used the 
4M Q ZAMS model, and assuming Q = 0.2 (black dots) and 
Q = 0.4 (red dots) we have calculated <j-modes in the frequency 
range between id = 0.1 and 1. 



Figure 12. Coupling coefficient e^j. between even \m\ = 1 g- 
modes associated with Aom and A4 m in panel (a) and between 
those associated with A2 m and A4 m in panel (b). Here, we have 
used the AMq ZAMS model, and for Q = 0.2 we have calculated 
g-modes in the frequency range between Gi = 0.1 and 1. 



cases with and without the effects of the centrifugal force on 
the modes, and a slight difference in the derivatives may 
cause a large difference in the locations of g-mode crossing 
between the two cases having similar gradients duin/dQ.. 

Fig. ll5l shows a mode crossing between two g-modes as- 
sociated with \2m and \i m for m — 1 which has a broad 
interaction range in £l. The broad interaction is possible if 
the oscillation frequencies of the two modes have similar de- 
pendence on the rotation frequency (i.e., similar dun/dfl) 
around the crossing, and the coupling coefficient is suffi- 
ciently large. The latter requirement is obviously met for 
the crossings between A2m and A4 m (?-modes (see panel (b) 
of Fig. [12]). 



The frequency dependence on the rotation frequency 
can be estimated from equation (|12[1 , Since A jm is propor- 
tional to v 2 for sufficiently large v (except for prograde sec- 
toral modes with Aom and m < 0), then, ur becomes pro- 
portional to y/2Q/n independent of j (or I). Therefore, most 
of the tJR-fi curves become closer to parallel to each other 
so that broad interactions such as the one shown in Fig. [15] 
would occur frequently as Q. increases. This explains broad 
and large peaks of uj\ for tesseral modes with I > \m\ in 
the panels (b) and (d) of Fig. [7] and for retrograde sectoral 
modes with I = m with m > in the panel (b) of Fig. [8] 

The prograde sectoral modes are exceptional. Since Aom 
with m < is nearly constant for v > 1, they experience 
only narrow interactions with tesseral modes. This explains 
why peaks for the prograde sectoral modes seen in the panel 
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1=1 and i = 3 calculated as solutions to the dispersion relation, 
eq. (fl3]>. for the 4M Q ZAMS model. To be compared with Fig. [9] 
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(a) of Fig. [8] are narrower than those in the other cases. The 
bumps seen for Q > 0.2 can be understood as swarms of 
numerous weak and narrow interactions. 



6 CONCLUSION 

We studied the stability of low degree g-modes in uniformly 
rotating main-sequence stars with masses of 4M and 5Mq , 
by taking into account the effects of both the Coriolis force 
and the centrifugal force, usin g the method of calculation 
given by iLee fc Baraffd (|l995l ). From the analysis treating 
f2 as a small parameter we found that a slow rotation has 
destabilizing (stabilizing) effect on high radial-order retro- 
grade (prograde) g-modes, although the effects for relatively 
low-order modes are very small or absent. This effect can be 
understood from the relation between period change due to 
rotation and the optimal period for the kappa-mechanism 
at the Fe-opacity bump. 

Calculating eigenfrequencies of low degree <?-modes at 



selected rotation frequencies, we found that, on the other 
hand, rapid rotation tends to stabilize high radial-order ret- 
rograde <?-modes, and the stabilizing effect appears stronger 
for less massive stars with lower effective temperatures. 

Obtaining the eigenfrequency of a g-mode of a degree 
£ (associated with a Xkm) as a function of the rotation fre- 
quency f2, we found that it experiences mode crossings with 
g-modes of different I's (associated with different A fe / m s). At 
a mode crossing between two modes with the same m and 
parity (even or odd) a coupling occurs because these two 
modes are not independent in a rotating star. If an unstable 
mode crosses a damped mode, the former can be stabilized 
around the crossing, or vice versa. At a large rotation fre- 
quency, low I g-modes, retrograde ones in particular, tend to 
be stabilized by mode couplings with larger £ modes which 
are strongly damped. The prograde sectoral modes (£ = — m 
and m < 0) in rapidly rotating stars receive exceptionally 
weak damping effects from mode couplings, because the £1 
dependence of their frequencies is very different from the 
other modes. 



Stability of g-modes in rotating B-type stars 13 



0.41 



1NJ 

m 



0.40 



0.39 



1 


1 

A 




// 




yy — 








yy 




yy 




yy 








yy 




^y 




y? 






AT 
















// 
















— 




</ 

1 





0.13 



0.14 



0.15 



2E-5 



£V(GM/R 3 ) 1/2 



1E-5 





1 

\ 

N 




1 




s 

\ 


S 

/ 












1 1 



5.13 



0.14 



0.15 



Q/(GM/R 3 ) 1/2 



Figure 15. Complex ui of even retrograde (m=l) g-modes with 
£ = 3 (solid line) and I = 5 (dashed line) calculated as solutions 
to the dispersion relation for the 4Mq ZAMS model. 



tained adiabatic g-mode spectrum of two-dimensional rotat- 
ing polytropes. We hope that it will become possible in the 
near future to calculate the frequency spectrum and stability 
for g-modes of two-dimensional evolutionary models. 
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We derived a dispersion relation useful to study mode 
coupling including the radiative damping effect under the 
quasi-adiabatic approximation. In the dispersion relation the 
strength of a mode coupling is controlled by a coupling co- 
efficient, tjk between two g-modes belonging to Xj m and 
Afe m (j > k); a larger ejk leads to a stronger coupling. The 
value of €jk is larger when k ^ and when 1/o;r and Q 
are larger, and it tends to be larger for retrograde modes 
than for prograde modes. By using the dispersion relation 
for various combinations of modes, we could reproduce all 
types of mode couplings. 

To include the effects of the rotational deformation of 
the equilibrium structure, we employed the Chandrasekhar- 
Milne (C-M) expansion, where the deformation is assumed 
to be proportional to f2 2 P2(cos 6). This approximation 
would not be very accurate for a star rotating at a nearly 
critical rate. One way to avoid th e problem would b e to 
employ two-dimensional models, as lBallot et al.l l|2010l ) ob- 



